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Abstract 

Using kaon secondary interactions in the material of the Belle detector, we search for both 
inclusive and exclusive production of the G(1540)"'". We set an upper limit of 2.5% at the 90% C.L. 
on the ratio of the 0(1540)^ to A(1520) inclusive production cross sections. We also search for 
the G(1540)^ as an intermediate resonance in the charge exchange reaction K~^n pl^s- 
upper limit of < 0.64MeV at the 90% C.L. at me+ = 1.539 MeV/c^ is set. These results are 
obtained from a 397 fb~^ data sample collected with the Belle detector near the T{4S) resonance, 
at the KEKB asymmetric energy e'^e~ collider. 

PACS numbers: 13.75.Jz, 14.20.Jn, 14.80.-j 
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I. INTRODUCTION 



The observation of the 0(1540)"'^ pentaquark, an exotic bound state with the quark con- 
tent uudds, is one of the most puzzhng mysteries of recent years (see [1] for an experimental 
overview). The ©(1540)"'" was first observed in exclusive reactions at low energy [1]. Later 
several groups reported observation in inclusive reactions at higher energies [1]. Conversely, 
other experiments at high energies do not see the 0(1540)+ pentaquark although they do 
observe significantly larger yields of conventional hyperons than seen at the experiments 
that observe the 0(1540)+. In order to resolve this discrepancy it is frequently assumed 
that pentaquark production decreases rapidly with energy. A high statistics experiment at 
low energies is therefore important. In order to achieve this, Belle utilises the small fraction 
of tracks that interact with the material of the inner part of the detector. These secondary 
interactions are used to search for pentaquarks. Particles produced in e+e~ annihilation at 
Belle have quite low momenta; the most probable kaon momentum is only 0.6 GeV/c. 

Results from two analyses are presented. In the first, we search for inclusive production of 
the 0(1540)"^ via the KN 0(154O)"^X, 0(1540)"^ pK^ process, using the signal from 
inclusive A(1520) production as a reference. In the second, we search for exclusive 0(1540)+ 
production in the charge exchange reaction K'^n 0(1540)+ — > pKg. For this search, the 
yield of charge exchange reactions is used as a reference, allowing a direct comparison to the 
results of the DIANA experiment [2] . 

II. DETECTOR AND DATA SET 

These studies are performed using a 357 fb"^ data sample collected at the T(45') resonance 
and 40 fb"^ at an energy 60MeV below the resonance. The data were collected with the 
Belle detector [3] at the KEKB asymmetric energy e+e~ storage rings [4]. 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a sili- 
con vertex detector (SVD), a 50-layer cylindrical drift chamber (CDC), an array of aero- 
gel threshold Cherenkov counters (ACC), a barrel-hke array of time-of-flight scintillation 
counters (TOF), and an array of CsI(Tl) crystals (ECL) located inside a superconducting 
solenoidal coil that produces a 1.5T magnetic field. An iron fiux return located outside the 
coil is instrumented to detect muons and mesons (KLM). Two different inner detector 
configurations were used. For the first sample of 155 fb~^, a 2.0cm radius beampipe and a 
3-layer silicon vertex detector (SVDl) were used; for the second sample of 242 fb~^, a 1.5 cm 
radius beampipe, and a four- layer silicon vertex detector (SVD2), and a small-cell inner drift 
chamber were used [5]. 

A GEANT [6] based Monte-Carlo (MC) simulation is used to model the production of 
secondary pK pairs, and to determine the detector resolution and acceptance. 

III. SELECTION OF SECONDARY pK PAIRS 

The analyses are performed by identifying pK~ . pK^ and pKg produced at secondary 
vertices. Charged particle candidates are required to be positively identified based on the 
CDC {dE/dx), TOF and ACC information. In addition, proton and charged kaon candidates 
are removed if they are consistent with being electrons based on ECL, CDC and ACC 
information. Kg candidates are reconstructed from tt+tt" pairs that have masses within 
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ilOMeV/c^ of the nominal K^g 

mass (3cr window). Additional selection requirements are 
imposed on the quality of the Kg vertex, on the impact parameters of the daughter tracks 
and on the angle between the momentum and the direction from the interaction point (IP) 
to the vertex. 

The proton and kaon candidates are required to have an origin that is displaced from 
the IP. The pK~ , pK~^ and pK^ vertices are fitted and those with a radial distance 1 < 
i? < 11 cm are selected. Additional criteria on the quality of the pK vertex are applied. We 
consider secondary pK pairs only in the central part of the detector —0.74 < cos 6' < 0.9, 
where 6 is the polar angle of the secondary pK vertex. The distributions of the secondary 
pK^ vertices in the xy plane are shown in Fig. 1 for the SVDl and SVD2 data samples, 
where the 2;-axis passes through the IP and is antiparallel to the beam. The beam pipe. 
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FIG. 1: Distribution of reconstructed secondary pK'^ vertices in the Belle detector for the SVDl 
(left) and SVD2 (right) data samples. 

the SVD layers, the SVD cover and the CDC support cylinders are clearly visible. 

IV. SEARCH FOR INCLUSIVE 0(1540)+ PRODUCTION 

The mass spectra for pK~ and pK^ secondary vertices are shown in Fig. 2. We apply an 
additional selection requirement on the angle between the pK momentum and the direction 
from the IP to the pK vertex d(j) < 1 rad. In the pK~ sample we reject A —>■ p7r~ and Kg 
TT+TT" decays misidentified as secondary pK~ vertices. No significant structures are observed 
in the rUpj^o spectrum, while in the rripK- spectrum a A(1520) signal is clearly visible. The 
pK~ mass spectrum is fitted to a sum of a A(1520) signal function and a threshold function. 
The signal function is a D-waye Breit-Wigner shape convolved with a detector resolution 
function, determined from MC. The detector resolution function is parametrized by a double 
Gaussian with widths of 2 MeV and 5 MeV with approximately equal contributions from each 
Gaussian component. The A(1520) parameters obtained from the fit are consistent with the 
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FIG. 2: (a) Mass spectra of pK (points with error bars) and pKg (histogram) secondary pairs. 
The fit is described in the text, (b) Momentum spectrum of the A(1520). 



PDG values [7]. The A(1520) yield, defined as the fit function signal component integrated 
over the 1.48-1.56 GeV/c^ mass interval (2.51), is (4.02 ± 0.08) ■ 10^ events. 

The pKg mass spectrum is fitted to a sum of a 0(1540)"'" signal component and a third 
order polynomial. The ©(1540)"'" is assumed to be narrow and its shape is determined by 
the detector resolution function, which is again a double Gaussian with similar parameters 
to those used for the pK~ mode. The MC resolution is checked against data using the 
fl~ —>■ AK~ signal, which has a topology similar to the secondary pKg pairs. It is found 
that the Q width is 4% percent larger in data than in MC, and this correction factor is 
applied to the mass resolution of secondary pKg vertices. For m — 1540MeV/c^ the fit 
yield is 58 ± 129 events. Using the Feldman-Cousins method of upper limit evaluation [8], 
we find N < 270 events at the 90% C.L. The upper limit is below 320 events for a wide range 
of possible G(1540)+ masses. We set an upper limit on the ratio of G(1540)+ to A(1520) 
production cross sections: 



iV, 



6(1540)+ epK- g(A(1520) pK- 

iVA(i520) epKO i3(e(1540)+ ^ pKl)B{Kl 



7r~^7r ) 



— < 2.5% at the 90% C.L. 



pK^s) = 25%. We take B(A(1520) ^ pR-) = 
[7]. The ratio of efficiencies for 0(1540)"'" pKg 



It is assumed that ,6(6(1540)+ - 
ii3(A(1520) ^ NK) = i(45 ± 1) 
and A(1520) — > pK~ of about 41% is obtained from MC simulation assuming that the two 
processes have similar kinematics. Our limit is much lower than the results reported by 
many experiments that observe the 0(1540)+. For example it is two orders of magnitude 
lower than the value reported by the HERMES Collaboration [9]. 

The momentum spectrum of the produced A(1520) is shown in Fig. 2 (b). This spectrum 
is obtained by fitting nipK- in momentum bins and correcting for the efficiency obtained from 
the MC. If a A(1520) is produced as an intermediate resonance in the elastic scattering of a 
K~ on a free proton, then its momentum is around 400MeV/c. As such, the observed hard 
momentum spectrum of the A(1520) confirms that they are produced in inelastic interactions. 
We find that non-strange particles do not produce A(1520) since the secondary pK~ pairs 
are accompanied by an additional K'^ from the same vertex in only about 0.5% of events. 
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The projectiles that can produce A(1520) are K~ , K^, Kl and A. The A needs a momentum 
of about 1.8 GeV/c to produce the A(1520). The number of such A's in e+e" annihilations 
is small. The observed number of A(1520) can not be produced by A projectiles even if 
the total KN cross section is assumed to be saturated by inclusive A(1520) production. We 
therefore conclude that our A(1520) signal is due to kaon interactions, with a contribution 
of A interactions no larger than a few percent. 

V. SEARCH FOR EXCLUSIVE 0(1540)+ PRODUCTION 

Possible exclusive pentaquark production is studied using the K^n — > pKg reaction, 
searching for the ©(1540)+ as an intermediate resonance. Since the projectile is not re- 
constructed, it is not possible to distinguish this reaction from the reactions K^p pKg, 
K^p pK^s- ^^^^ inelastic reactions with a ti^ or undetected tracks in the final state. There- 
fore we apply selection criteria which suppress the contribution of the inelastic reactions in 
the sample of secondary pKg pairs (see section V E) , and determine the contribution of the 
charge exchange reaction indirectly, as described below in this section. The cross section 
of the K^p 0(1540)"'" pKg reaction is expected to be similar to the cross section of 
the K'^n 0(1540)+ pKg reaction [10], however the fiux of Kg is decreased by decays 
in flight and we conservatively neglect the contribution of this reaction. The K^lp pKg 
reaction can produce a dip in the 0(1540)+ region because of the interference of resonant 
and nonresonant amplitudes. However, the size of the dip is expected to be negligible [10]. 

The number of charge exchange reactions can be estimated in a straightforward way 
using the known flux of primary i^+, the reaction cross section, cr^^, the amount of 

material, M, and the reconstruction efficiency for the secondary pK^ pair, Cp^o, and taking 
into account nuclear effects: 

N'^^inipK) = J (J'^^ M €pKO B S{En, \pf\) S{^/^-mpK) P dE^ d^Pr dpK+ dR d9, (1) 

where B is the product of branching fractions B = B{K^ — > Kg) B{Kg — > 7t~^t:~), 
S{E]\f, \pf\) is a nuclear spectral function which is a joint probability to find in a nucleus 
a nucleon with energy E^- and Fermi momentum pp, s = {Ek+ + -E'at)^ — {pk+ + PfY is 
the centre of mass (cm.) energy of the reaction squared, Ek+ is the energy and Pk+ is the 
momentum of the projectile, rUpK is the mass of the produced pair, P is the probability that 
the produced pair is not rescattered in the nucleus and R and 9 are the radial distance and 
polar angle of the secondary vertex, respectively. 

However, it is difficult to accurately estimate the systematic errors in this calculation 
because M and e^^o are complicated functions of the coordinates and the estimation of 
S and P is model dependent. This problem is solved by reconstructing the decay chain 
D*^ — > £)°7r~, — > X+TT" for events where a K'^ interacts elastically in the detector 
material. The reconstruction procedure is explained in detail in section VA. The yield of 
such decays, N^* , is expressed as 

N^D*{mpK) = J ^D* t^'' M epK+ S{En, |p>|) ^{^/^ - rripK) P dE^ d^Pr dpK+ dR d9, (2) 

where is the fiux of K^ originating from the selected D*~ decay, cr*^' is the cross 

section for elastic X+p — > pK'^ scattering and epK+ is the efficiency of reconstructing the 
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secondary pK^ vertex. The nuclear suppression P is assumed to be the same for pK'^ 
pairs and pK^ pairs produced in the charge exchange reaction. We express A^'^'^ in terms 
of A^lJ*. The expressions are simphfied by the facts that the products Me^^o and MepK+ 
are approximately independent of 6 for the central part of the detector; the ratio ^pK%l^vK^ 
is approximately independent of the secondary pair momentum, and the nuclear 

suppression P is approximately independent of \P'pk\ (see section VA). We approximate 
the nuclear spectral function as S{Eiy, \pf\) = ^{\Pf\) ^(-E'Ar — fUppD), where the function 
is defined so that = corresponds to the maximum of S{Ej^, \pf\)- We 

obtain 



where 

^^^jmpx) _ J ^^^{\pk+\) W{\pf\) 5{y/s-mpK) epK+jnipK, \Ppk\) d^pF dpK+ 
^^tirUpK) ^ ^^t{\pK+\) W{\pf\) S{Vs-mpK) epK+{mpK, \Ppk\) d^pp dpK+ 



(4) 



The p.d.f. of the Fermi momentum distribution l^(|i7p|) is determined from data in sec- 
tion VA. The form of the is discussed in the same section. The integrations are 
described in section VB. Equation (3) provides the basic formula to determine the yield of 
the charge exchange reaction in our data sample. The formula (3) would be even simpler if 
we could use projectiles from D*^ decays that undergo charge exchange in the detector 
material. However, the fraction of such events is very low because of the relatively small 
cross section, and the low reconstruction efficiency and branching fraction. 



A. Determination of iV|}, 

To determine the number of D*^ D^tt^, — > K'^n~ decays for which a in- 
teracts elastically in the detector material, the four-momentum of the interacting is 
reconstructed based on the information available from the produced secondary pK^ pair. 

If a secondary pK pair is produced in a quasi-elastic reaction, four-momentum conserva- 
tion implies 

Er + En = EpK, (5) 
Pk+Pf=PpK- (6) 

The nucleon energy, Ej^ is approximated by [11] 

Er,^mr,-2e-^, (7) 

where m^r is the nucleon mass and e ~ 7 MeV is the nucleon binding energy. This approx- 
imation is valid for the high \pf\ part of the nuclear spectral function. Another possible 
approximation is Ei\f = mj^ — E^i, where Eji ~ 26 MeV is the average removal energy of the 
boimd nucleon. We find that both approximations give a similar resolution in the mass 
(see below) and our result is independent of the choice. The quantities EpK and ppK are 
measured; taking into account the primary and secondary vertex constraints, Eqs (5)-(7) 
can then be solved iteratively. 
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The iterative process is started from Eq. (7) where some average value of \'Pf\ is substi- 
tuted. Then the energy of the projectile Ek is determined from Eq. (5). In the next step 
the projectile momentum is determined from its absolute value = yj E]^ — t^K, and 
the flight direction obtained from the primary and secondary vertex constraints, taking into 
account the bending of the track in the magnetic fleld. The value of \pf\ is then determined 
from Eq. (6) and the iteration loop is closed by substituting the obtained \pf\ into Eq. (7). 

The projectile four-momentum is determined for all secondary pK^ pairs. The resulting 
projectile candidates are then combined with all the 7r~ candidates in the event to form 

candidates; the D'^ candidates are combined with all the remaining tt^ candidates to form 
D*'^ candidates. The tt" candidates are required to be positively identified based on the 
CDC (dE/dx), TOP and ACC information and to originate from the vicinity of the IP. We 
reject vertices with additional tracks and require 50 < < 300MeV/c. The lower bound 
on \pp\ is used to reject interactions on hydrogen (see below). Events in a ±2MeV/c^ {3a) 
window in ArriD* = uid* — rn-^o are selected and the mass of the daughter D° candidates is 
plotted (see Fig. 3 (a)). A signal of 470±26 D° events with the mass consistent with the PDG 




FIG. 3: (a) Invariant mass of the selected candidates with the K'^ reconstructed via a secondary 
pK~^ vertex, (b) Sideband subtracted distribution of \pf\ for the selected candidates, fitted to 
the oscillator model expectations. 

value and a mass resolution of 16MeV/c^ is observed. The D° signal and sideband regions 
are selected as \mK+n- ~ '>t">-do\ < 50MeV/c^ and 60 < IrriK+n- ^ fT^Dol < llOMeV/c^, 
respectively. The sideband subtracted distribution of \pf\ is shown in Pig. 3 (b). The 
peak near zero is attributed to interactions on hydrogen, which is present in the detector 
material. If En = ttin is used instead of Eq. (7), this peak is found at zero, as expected for 
interactions with a free proton. The \pf\ spectrum is fitted to the parametrization, expected 
in the oscillator model [12]: + 4/3(|Pf|/Po)^] exp(— |pVP/pq). The value for the model 

parameter returned by the fit, po — 115 ± 4MeV/c, is comparable to those obtained from 
other measurements of \pf\ distributions [12]. 

We find that the fraction of D*~ events in the hydrogen peak is roughly independent 
of pair momentum. Since this fraction is inversely proportional to P (there can be no 
rescattering in hydrogen), we conclude that P is also roughly independent of pair momentum, 
which is important to perform the integration in Eq. (4). 
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The m-jjo spectra are fitted in nipK bins to determine the D*~ yield. The fit function is 
comprised of the sum of a Gaussian and a first order polynomial The Gaussian mean and 
width and the polynomial parameters are all fioated in the fit. The number of D*~ mesons 
in the 6(1540)+ region is 24 ± 7 per SOMeV/c^ bin. 

B. Determination of $^^/$^*^ 

The fiuxes of primary K~^^s and K~^^s from D*~ are determined from data. Primary 
candidates are required to originate from the vicinity of the IP and to be positively identified 
based on the CDC (dE/dx), TOP and ACC information. D*^ candidates are selected by 
combining the kaon candidate with pion candidates, in the same way as described above. 
Correction for reconstruction efficiency and contamination from other particle species is 
performed using MC that is calibrated from data. The integration in Eq. 4 is performed 
using a Monte Carlo technique. The nucleon Permi momentum pp is assumed to be isotropic 
relative to the projectile momentum pk- The fact that the secondary pK pair efficiency 
depends on the momentum results in a correction of about 3% in the ratio. The flux ratio at 
f^pK = 1-539 GeV/c^ is equal to 850 ± 20; the uncertainty is dominated by the assumption 
concerning the relation between Ej^ and \pf\- For this measurement we use about 20% of 
the data sample distributed uniformly over the running period. 

C. Determination of a'^^/a^^ 

The cross sections cr^^ and a^^ are obtained from published data [7, 13]. The data are 
fitted with polynomials in rUpK, and the ratio of the fitted functions is used to obtain a'^^/a'^^. 
The value of the ratio is 0.35 ± 0.02 at rUpx = 1.539 GeV/c^ and rises with rUpK- Errors are 
assigned based on the typical experimental errors in the region of interest. 

D. Determination of ^pK^/^pK+ 

Monte Carlo simulations are used to estimate the ratio epKO^/epK+- The angular distribu- 
tion in the reaction cm. frame is assumed to be uniform as expected for low energy elastic 
K^p scattering and for 0(1540)+ pKg decay. We consider the following sources of sys- 
tematic uncertainty: Kg, and secondary pi^T reconstruction efficiency (7%), uncertainty 
in material description (5%), uncertainty in the description of the reaction kinematics (5%), 
and the MC statistical uncertainty (5%). The ratio of efficiencies at rUpK — 1.539 GeV/c^ is 
(43 ± 5)%. Here all the uncertainties are added in quadrature. 

E. Upper limit on G(1540)+ yield in exclusive reaction 

To suppress the contribution of inelastic reactions in the sample of secondary pKg pairs 
we reject vertices with additional tracks and require 50 < \pf\ < 300MeV/c. The lower 
bound on \pp\ is used to reject interactions on hydrogen, which do not contribute to the 
charge exchange reaction. The effect of angular cuts used by DIANA to suppress rescattering 
is checked, and found not to suppress background significantly when the \pf\ cut is applied, 
and not to improve the sensitivity. The pK^ mass spectrum and the expected yield from the 
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charge exchange reaction is shown in Fig. 4 (a). 



The statistical and systematic uncertainties 
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FIG. 4: (a) Invariant mass spectrum for secondary pKg pairs (small dots with error bars) and 
expected yield of the charge exchange reaction per 2MeV/c^ (open dots). A fit to a third order 
polynomial is shown with a dashed curve. The 0(1540)+ contribution expected from the DIANA 
result is shown with solid line, (b) The yield of 0(1540)+ from the fit, expressed in terms of the 
resonance width. The open dots correspond to the upper limit at the 90% C.L., obtained using 
the Feldman-Cousins method. The square with error bars indicates the current PDG value for the 
G(1540)+ width. 

on A^'^^ are added in quadrature. In the 6(1540)+ region we expect (1.03 ± 0.36) ■ 10^ 
charge exchange events per 50MeV/c^ bin. The mp^o distribution is fitted to a third 
order polynomial and a signal p.d.f. positioned at various values of nipj^o. The fit finds 
A^e+ = -11 ± 59 candidates at = 1.539 GeV/c^. The ratio of the 6(1540)+ yield to 
the charge exchange reaction yield can be expressed in terms of the 6(1540)+ width (see for 
example [14]): 



F« 



a 



ch 



e+ 



Am, 



iV^h 107 mb Bi Bf 

where Bi and Bf are 6(1540)+ branching fractions into initial and final states, B^ = Bf = 
0.5, and Am is the mass interval of pKg pairs used to determine A^'^^. The resulting values of 
Fe+ are shown as a function of mpx in Fig. 4 (b). Also shown in Fig. 4 (b) are the 90% CL. 
upper limits, obtained with the Feldman-Cousins method [8], and the current PDG value 
for the 6(1540)+ width [7], which is based on the DIANA result. A similar width has been 
inferred from a reanalysis of K~^d scattering [15]. It is assumed that nuclear suppression for 
the 6(1540)+ is the same as for nonresonant pKg pairs. 



VI. CONCLUSIONS 

Using kaon interactions in the material of the Belle detector, we searched for both inclusive 
and exclusive production of the 6(1540)+. No 6(1540)+ signal was found in the sample of 
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secondary pKg pairs. For inclusive production we set the following upper limit: 



a(KN e(1540)+X) ^ , ^ ^ 

For exclusive production we find 

r(K+n e(1540)+ ^ pK^s) < 0.64 MeV at the 90% C.L. 

at me+ = 1.539 MeV/c^. This upper limit is below the current PDG value of 
r = 0.9 ±0.3 MeV, and below 1.0 MeV for a wide interval of possible 0(1540)+ masses. 
This measurement uses a sample of low energy kaon interactions and allows for a direct 
comparison with the DIANA result. With similar sensitivity, our results do not support 
their evidence for the 0(1540)"'". 
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